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Abstract 

We consider lattice implementation of the recently proposed gauge invariant defi- 
nition of the monopole charge. Because of the lattice discretization the algorithm 
gives rise to specific lattice artifacts and an effective Ising model. The Ising-model 
problem might in principle be solved and we discuss the role of the Maximal Abelian 
gauge in this respect. The lattice artifacts are much more difficult to deal with since 
they are mixed up with monopoles thus obscuring the physical observables. Never- 
theless, it is possible to extract the density of physical monopoles which seems to 
scale correctly towards the continuum limit. 



Introduction 



In this paper we implement the gauge invariant monopole construction M on the lattice. 
Since we use heavily the results of Ref. [Ij let us recall first the essentials of the monopole 
definition in the continuum limit. 

One starts with consideration of the fundamental Wilson loop W(T), ^TrW(T) = 
cos (p(T) calculated on a closed contour C(t), t e [0;T]. For a given W(t) it is always 
possible to find a family of spin coherent states | n(t) ) such that 

e^ \n{t)) = W{t)\n{0)), n(T) = n(0) . (1) 

Eq. (|TJ) naturally assigns the unique (up to the sign, see below) vector n(t) to each point on 
the contour C. In terms of n{t) the phase angle <p{T) of the Wilson loop can be expressed 
as 

(p(T) = - An + - n-dnx dn , (2) 
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where Sc is an arbitrary surface spanned on C and A denotes the tangential to C com- 
ponents of the gauge potentials, A a (t) = A^(x(t)) x^{t). As is clear from Eq. ([]]), the 
initial state | n(0) ) is an eigenstate of the full Wilson loop. Since W(T) has two eigen- 
states | ± n(0) ) ('spin up' and 'spin down'), with the corresponding families ±n(t) and 
phases ±<p(T), there is a sign ambiguity in Eq. ([!]). 

One considers next an arbitrary smooth closed surface Sp hys and covers it with a set 
of infinitesimal patches of area Sa x . The phase evaluated on each patch is: 

5ip x = - { d A (nA) + n ■ dn x dn} 5a x . (3) 

One can show that the usual assumption of the gauge fields continuity implies existence of 
a smooth field n x e S^ h . Therefore it is legitimate to integrate Eq. (^) over the surface 

C2 . 

q = — <b 8<p — — <f> {d A (nA) + n ■ dn x dn} 5a , (4) 
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and get a gauge invariant definition of the monopole charge contained inside Sp hys . Note 
that for a smooth S^ hys and for continuous gauge fields the sign ambiguity mentioned 
above is not important. In fact, it reduces to the freedom to globally change the sign of 
n G Sp hys which in turn is equivalent to changing the sign of q, see Eq. ([|). 

One can show [p]] that at the classical level it is the singular Wu-Yang monopole 
solution which corresponds to a non-trivial q. In other words, Eq. (|J) provides us with 
a gauge invariant topological definition of the Wu-Yang monopole which can be used at the 
quantum level as well. Moreover, the quantum generalization of the Wu-Yang monopole 
may look quite different from its classical counterpart. Indeed, the classical Wu-Yang 
solution is unstable 0. Thus, the semi-classical approximation would be inconsistent. 
However, the instability of the Wu-Yang solution implies only that monopoles (|J) are 
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irrelevant at short distances r <C Aqcd- As far as the infrared region is concerned the 
dynamics of the monopoles (f|) cannot be treated analytically (see also Ref. ||). One 
concludes, therefore, that at the quantum level the field configurations with a non-trivial 
charge q have little to do with the classical Wu-Yang solution both at short and large 
distances. Moreover, there seem to be no analytical tools to investigate their dynamics. 

Below we discuss the monopole definition on the lattice. Lattice implementation 
of the continuum ideas is not always straightforward. A particular example is provided 
by instantons: while in the continuum limit instantons are relatively simple objects, their 
lattice definition is technically much more involved 0. Still, in principle, one may try 
to interpolate the lattice gauge fields and apply the continuum ideas to the interpolated 
gauge potentials. A drawback of this approach is that the interpolation is always not 
unique. 

Because of the lattice discretization one cannot completely avoid interpolations while 
implementing Eq. ([|) on the lattice. However, our interpolation procedure is unique 
and not applied to the gauge fields themselves. Instead we interpolate the states | n ) in 
between neighboring lattice cells thus getting a continuous field n as required by Eq. (f|). 
This enables us to construct integer valued conserved monopole currents on the lattice. 

There are two disadvantages of our approach both of which we discuss in detail later: 

i) The sign ambiguity mentioned above which is harmless in the continuum reemerges 
on the lattice in the form of an effective Ising model. In turn, the problem of finding the 
gauge invariant monopoles on the lattice becomes equivalent to finding the global minima 
of the corresponding Ising action, which is non-trivial because the model is frustrated. 
What is surprising, however, is that in the Maximal Abelian gauge the ground state of 
the Ising model is approximated very well, although we do not fix any particular gauge 
neither in our construction nor in the numerical simulations. 

ii) The interpolation of spin states in between neighboring lattice cells enforces a 
change of the lattice geometry. Namely, one has to consider additional three dimensional 
cells on the lattice which can contain gauge invariant monopole charges. We argue below 
that these charges are pure lattice artifacts which disappear in the limit of infinitesimal 
lattice spacing. The artifact 's contribution to the physical observables must be carefully 
subtracted which is a non-trivial problem because of complicated lattice geometry and 
monopole-artifact mixing. Nevertheless, for simple quantities like monopole density the 
monopole-artifact separation turns possible. When this is done, the density of physical 
monopoles scales correctly towards the continuum limit. 

1. 'Stokes Theorem' 

Complexity of the lattice implementation of Eq. (|]) might already be visualized from 
that equation itself. Indeed, the effective 'Higgs field' n is unusual since it is defined not 
on space-time points (lattice sites) but on elementary two-dimensional cells (plaquettes) . 
Moreover, the field n is defined [|] in terms of the field strength tensor F® u which itself is 
badly defined on the lattice. Indeed, in the continuum limit the Wilson loop calculated on 
a plaquette reduces to a particular component of F^ v regardless of the initial point from 
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which we started evaluation of Wilson loop. But for a finite lattice spacing the Wilson 
loop depends drastically on the starting point. Because of this, the spin state which we 
would like to assign to a given plaquette also depends on the plaquette definition. 




Fig. 1 

The situation becomes even worse if we consider a common link of two plaquettes. 
The spin states are to be constructed on each plaquette separately and they have nothing 
in common, in contrast with the continuum case. To get insight into the problem, let us 
consider an arbitrary contour Cq and dissect it in points A and B into two smaller pieces 
Ci and C 2 (Fig. 1). The spin states are introduced for each contour via Eq. ([]]) taking for 
definiteness the 'spin up' initial states. In this way one gets three different families Hi(t) 
associated with each contour Cj. Suppose that we know the phases <£>j of the Wilson loops 
calculated on Ci. What is the relation between them? 

The answer may be found by considering the image of each contour under the map: 
H i : Ci — > Cl° loT . Note that the two triples {n^A)}, {n^B)}, i = 0, 1, 2 are distinguishable 
since the relative position of the vectors in each triple is gauge invariant. Moreover, there 
exists a simple interpretation of the original Wilson loops in terms of C^ olor . Namely, Wi(t) 
describes periodic motion of a point-like charged particle on S^ olor since every closed path 
(jcoior j g aS s OC i a ted with a unique phase factor. Therefore, in terms of Cf olor the problem 
becomes essentially Abelian and one could have concluded that (fo is equal to the sum 
of phases calculated on the pre- images of the closed paths I, II, VII (see Fig. 1). But 
there is a loophole in this argumentation: we do not know the SU(2) matrices gij which 
connect the points n» — > Hf 

g i3 {A) : Hi(A) -> H 3 {A) , 9ij {B) : n^B) -> H 3 (B) (5) 

and which are denoted by dashed lines on the Fig. 1. Clearly, g i3 - might be an arbitrary 
SU(2) elements provided that Eq. (|) holds true. 

Among various g G SU{2) which connect two given points ni,n 2 G S^, olor , g\fi\) ~ 
| n2 ) there is a special one. Namely, if g describes the motion n\ n 2 along the shortest 
geodesic connecting Hi, n 2 then g is uniquely defined: 

g = (nin 2 ) + ia a [n 1 xn 2 ] a - (6) 
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We call it the geodesic matrix below. Moreover, for the geodesic matrices the relation 
between the phases (fii becomes extremely simple. Indeed, one can readily verify that in 
this case the phase angles associated with cells V, VI, VII are vanishing while for the 
the cells III, IV they are equal to the oriented areas of the spherical triangles {fii(A)} 
and {raj (£?)}, respectively. Thus, we get the following simple equation which relates the 
phases of the Wilson loops: 

cp = ^ + cp 2 + 7 (A) + 7 (5), (7) 

7(iV) = oriented solid angle between no(N), Hi(N), ^(N) . 
Note that we always implicitly assume mod 2ir operation [1]. 



2. Single Plaquette Construction. 

We start the construction of monopoles @ on the lattice following the same strategy as in 
the continuum limit. Namely, the simplest object which we consider first is an elementary 
plaquette, see Fig. 2. 




a . 5 

g(A) 



n/A) 



n,(B) 



g(B) 



n 2 (A) 



n 2 (B) 



Fig. 2 



P Y 
Fig. 3. 



As is noted above for any finite lattice spacing the matrix U p = U1U2U3U4 is to be 
treated as a large Wilson loop. In particular, it is impossible to assign a single well defined 
spin state to a given plaquette. On the other hand, the Eq. (ffl) is still valid and in fact 
directly applicable. Namely, one can calculate the eigenvector of the matrix U p : 

(n \U p = e^(n | (8) 

and then propagate the initial state ( n Q | along the plaquette's edges: 

(n k \U k = e l<Pk (n k+ i\, k = 0, ...,3, (n 4 | = (n |. (9) 

In this way one gets four states ( n k \ sitting at the plaquette's corners. Moreover, each 
( fik I is an eigenstate of the plaquette matrix calculated starting from point k, e.g. 
(fii \U1U2U3U0 ~ (ni I etc. Therefore, the vectors are in a sense always parallel to 
the color magnetic field piercing the plaquette considered. Moreover, Eqs. (||,[]) imply the 
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following simple relation between the plaquette angle (p p and phases (fi associated with 
the corresponding links: 

3 

V P = ^Vk- (10) 

Note that we still have a freedom to choose an arbitrary sign of the initial state (H), 
the same as in the continuum limit. It is clear that the change n Q — > — n would also 
flip the signs of all other states: — > — n&. Thus, it makes sense to speak about the 
sign associated with each plaquette. While only a single plaquette is considered the sign 
ambiguity is in fact irresolvable. Note also that for pure diagonal link matrices all vectors 
n k collapse to a single point on the S^ olor . 

Next, consider a common link of two adjacent plaquettes, see Fig. 3. Since the spin 
states are constructed separately on each plaquette the states fii(A), n 2 (A) and r?i(-B), 
n 2 (B) are quite different in general. It is important, however, that the products n\(A) ■ 
712(A) and fii(B) ■ n 2 (B) are gauge invariant. Moreover, in case of diagonal lattice gauge 
fields (or for a gauge copy of the diagonal configuration) we have 

n x {A)-n 2 {A) = n^B) ■ n 2 (B) = ±1, (11) 

depending on the relative sign of initial states chosen on each plaquette. Therefore, there 
is a natural way to fix the relative sign on Pi and P 2 by the requirement that across 
the common link of two plaquettes the spins should be as smooth as possible. In fact, 
this approach is essentially the same as in the continuum limit. The only difference is 
that in the continuum we indeed do have a continuous distribution of spins while on the 
lattice the true smoothness of n x is not possible. The only exception is a particular case 
|ni(^4) • n2(^4)| = 1 for which everything is in fact trivial. 

Thus, on the lattice we inevitably have to interpolate the states n\(A) and n 2 (A) by 
connecting them with an appropriate SU(2) matrix g such that ( Hi \g ~ ( n 2 |- As we have 
argued in the previous section there is a natural unique choice of g, namely the geodesic 
interpolation, Eq. (|6]). This choice is singled out, in particular, since only in this case we 
do not introduce additional fluxes in the intersection of two plaquettes: the Wilson loop 
calculated on 2-cell (a(3^(5) is exactly unity no matter what the lattice gauge fields are. 

3. Monopole Charge on 3D Cube. 

The above construction allows to formulate the notion of the monopole charge for a single 
three-dimensional lattice cube. Namely, one has to apply the procedure of the previous 
section to all six plaquettes of a given cube, Fig. 4. 

Let us briefly summarize the essential steps of the construction. 

i) On each face a = 0, ...,5 of the cube one calculates the corresponding plaquette 
matrix U a and then constructs four vectors nj iCr , i = 0, 3 via Eq. (|9]) taking the initial 
state ra ,a as an eigenstate of U a with arbitrary sign. 
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ii) In order to get the smoothest distribution of n i>a on the cube boundary, one con- 
siders the following Ising-like action: 

S = - V, ZgZp n i>a ■ n jtj3 , (12) 

a/3 i,j 

where dynamical variables are the signs z a = ±1 associated with each plaquette, the 
first sum is over all neighboring 2-cells in the cube's boundary and the summation in i, 
j indexes extends over all neighboring spins. The action ( |T2"D is minimized with respect 
to all signs z a and, finally, for negative z a the corresponding vectors n ij0l are inverted 
— ► —n^a- It is clear that the ground state of the Ising model flT2"| ) corresponds to the 
smoothest possible distribution of ni ta . Moreover, for a single cube the ground state may 
be found exactly. 

in) Finally, the three vectors around each vertex of the cube are geodesically interpo- 
lated (dashed lines on Fig. 5). Note that the relative position of vectors in each triple, say 
{ni,p ^3,7 ^0,5) is gauge invariant and therefore the geodesic interpolation is a well defined 
prescription. 

The magnetic charge contained in a given cube is proportional to the net outgoing 
magnetic flux, cf. Eq. (f|). Namely, one sums up the phases ipi associated via Eq. (@) with 
all two-dimensional cells on the boundary of the 'extended' cube (Fig. 5, right). Note 
that due to the geodesic interpolation not all the 2-cells must be taken into account. For 
example, the flux emanating through the 2-cell {ni,/3, ^o,/3, ^0,7, ^3,7) is exactly zero. The 
final expression for the magnetic charge is therefore 

(13) 

where (p a is the plaquette angle (H) and (p v is the flux corresponding to the triangle (dashed 
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lines) at vertex v which in turn is equal to the oriented solid angle between vectors ni )Ql 
at the triangle's corners. 

Eq. ( |13|) has in fact a transparent meaning. The first term counts the flux emanating 
through the original lattice plaquettes and therefore it is the physical flux of the monopole. 
The second term is in a sense a correction due to the lattice coarseness. Indeed, there were 
two important requirements in the corresponding continuum construction : S% hys should 
be smooth and 5ip x is to be calculated on an infinitesimal contour 8C X . In fact none of them 
is realized on the lattice since the lattice cube for sure cannot be considered as a smooth 
surface and the lattice plaquette is not infinitesimal in any sense. In fact the second term 
in Eq. (|13|) just accounts for these facts by plugging the holes between the neighboring 
vectors n x . Let us emphasize that the appearance of the additional contribution, Eq. fllcf) , 
is entirely due to the lattice discretization. There are no similar term in the continuum 
expression (f|). 

However, there exists a particular case when the second contribution, Eq. (|13|), is 
zero. This happens when all the neighboring spins n^ a are aligned with each other. From 
Eq. ([11]) we conclude then that this corresponds to a pure Abelian configuration. Namely, 
there exists a particular gauge in which all the lattice gauge fields are diagonal. Moreover, 
Eq. ( |13D is then identical to the well known DeGrand-Toussaint definition of the Abelian 
monopole |J. 



4. Effective Ising Model and Lattice Artifacts. 

In this section we discuss some details of the above construction which are in fact crucial 
in the actual numerical implementations. For simplicity, we consider in this section three 
dimensional lattice only. Let us start from the Eq. (|T2| ) which encodes the way how to 
get the smoothest possible distribution of spins fij )Q .. In fact, it is quite important to find 
the global minimum of the action (0). Indeed, let us consider an almost Abelian gauge 
fields configuration for which 

^ w ± S a3 VijQ ,_ (M) 

Suppose that instead of the ground state of the model ([12]) we are considering a local 
minimum of the action (|12). After performing the corresponding inversions ni )Ql —>■ — ni jQ; 



for all the negative z a we get a configuration of vectors nj jCf almost all of which satisfy 
ni,a ~ ^ a3 - But since the minimum is only local there are some negative spins as well. 
Evidently, Eq. (|T^) when applied for such {n ita } still produces a non zero charge which 
however has nothing to do with reality since the contribution comes mostly from the 
second term of Eq. fljjD . 

Therefore, it is indeed important to find the true ground state of the model (O) 
before applying Eq. (|13|) . One expects, however, that in the real simulations it is almost 
impossible to find the global minima since for non-Abelian gauge fields the model flT2"|) is 
frustrated. Thus the concrete algorithm used to minimize flT2|) might be crucial. 

As a matter of fact, the naive local minimization of (0) is doomed to produce unphys- 
ical results for the monopoles fllTf). One may try to improve the iterative minimization 
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by preconditioning it with a particular gauge fixing of the original lattice gauge fields. 
In our numerical tests (see below) we have considered three particular gauge fixing pre- 
scriptions: minimal and center Landau and Maximal Abelian gauges which are defined 
via maximization of the following functionals 

£^ Trf/ " Ef^) 2 ' E lT*U l0 *Ufo>, (15) 
i i ^ ' i 

respectively. Our finding is that the Maximal Abelian gauge is indeed greatly favoredQ 
Namely, when the lattice fields are first MAA gauge fixed then for each a it is sufficient 
to pick up the nearest to the north pole family n i)Cx . The resulting distribution {fii ja } 
almost exactly corresponds to the ground state of fll2"f ). For other gauges the result is less 
pronounced: although Landau gauge fixing improves the local minimization the resulting 
minimal value of ([12]) is considerably larger than in case of the Maximal Abelian gauge. 

Let us note that in our numerical simulations none of the gauges (|T5|) was used. 
Instead, we have applied a variant of the simulated annealing algorithm to directly 
minimize the action (0). As one may expect the latter approach turns out to be superior 
to any kind of local minimization. 

Let us turn now to another major problem of the lattice implementation of Eq. (§). 
Namely, one can see that Fig. 5 implies in fact a change of the lattice geometry since 
this construction being iterated for all 3-cubes produces 'extended' cubes instead and 
introduces new three dimensional cells on the lattice. In particular, in the intersection of 
four cubes one gets the geometrical figure depicted on the Fig. 6 and a new 3-cell at each 
site of the original lattice, see Fig. 7. 




Fig. 6 Fig. 7 



Consider first the Fig. 6. For the plotted 3-cell the faces (a^c^i), i = 1, 2 consist of the 
geodesic matrices while all others were in fact already discussed. Namely we have shown 

1 See the next section for details on numerics. 
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that fluxes piercing the cells (&1&2C2C1), (01,02,61,62) etc. are exactly zero. Moreover, 
since there is the same matrix U on each link (0102), (6162) etc. one concludes that the 
images of squares (dibiCidi), % — 1, 2 under the maps n(cij), n(6«) etc. are isomorphic. Due 
to the properties of the geodesic interpolation the fluxes piercing (ai&iCidi) and [a^c^d^] 
are exactly the same. 

The overall conclusion about Fig. 6 is as follows. The magnetic charge inside the plot- 
ted 3-cell is identical zero. Moreover, the flux which goes from this cell to the neighboring 
cubes is zero either. The whole flux propagates unmodified along the link of the original 
lattice from the face (ai&iCidx) to the face (02^202), see Fig. 6. 

Turn now to the 3-cell on the Fig. 7. Its properties are easy to derive. Namely, at each 
vertex there is a corresponding vector Hi, i — 0, 11. It is convenient to parameterize 
the vertices as (x; fi, d M ; v, d u ) = (2; v, d v \ fi, d^), where ft, v are the unit vectors along the 
lattice axes and d M , d v are the corresponding shift signs (see Fig. 7). Moreover, the pairs 
of vertices 

(x;fi, d^v, d u ), (x; fi,d^X,d x ) V>, v, X 
are connected by a link which carries an SU (2) matrix g such that 

g : n(x] fx, d^, v, d v ) -> n(x; fi, d^; X, d x ) 

corresponds to the geodesic motion on S^ olor . Note that the relative position of all 
n(x; fi, d^; v, d u ) is gauge invariant and therefore the construction of the 3-cell on the 
Fig. 7 is well defined. 

From the above properties one concludes that fluxes emanating from the faces of the 
presented 3-cell are equal to the oriented solid angles between the corresponding vectors 
situated at the vertices of that face. Then it follows immediately that the magnetic charge 
inside this 3-cell is gauge invariant and in general non-zero. Thus, we get another type 
of monopoles which however live in the sites of original lattice contrary to the monopoles 
in Eq. (|T3|), which belong to the dual lattice, as usual. In fact, there is another striking 
property of the monopoles on the Fig. 7. Namely, none of the fluxes which this monopole 
emanates enter the original Wilson action of the lattice gauge models. Therefore, naively 
this monopole costs no action at all. 

One can readily convince oneself that the actionless monopoles are of the same origin 
as the second term in Eq. (0). They appear because of the violation on the lattice of 
the two major requirements of continuum theory (see the discussion after Eq. (|T3"D). We 
conclude, therefore, that the actionless monopoles are intrinsic to the lattice formulation. 
They disappear only in the academic al limit of such small lattice spacing that one can 
indeed identify the lattice plaquette with a particular component of continuum F^ u . Of 
course, in this limit the above single cube construction is not applicable any longer and 
one has to consider the extended monopoles 0. 

Thus, the actionless monopoles are in fact lattice artifacts which exist only because 
of the lattice discretization. When any physical observable is measured on the lattice 
the contribution of this objects must be carefully subtracted - the task which is quite 
non-trivial by itself. Indeed, even if D = 3 and the lattice geometry is well understood 
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(see Figs. 5-7) there seems to be no unique way to extract, e.g., density of the physical 
monopoles. In three dimensions the monopoles are point-like and one may naively expect 
that each artifact should form a small dipole-like structure either with 'physical' monopole 
or with another artifact. But this expectation is in fact wrong since the flux emanating 
from each artifact is not suppressed by the action. Therefore, e.g., artifact-artifact chain 
might be arbitrarily prolongated. 

In four dimensions the situation is even worse because of the complicated lattice ge- 
ometry. The monopole currents are of course conserved, but the conservation takes place 
not on the hypercubical lattice. As a result the monopoles may freely propagate from 
3-cubes to the sites and back on the original lattice. Thus, discrimination against the 
artifacts becomes ambiguous. Nevertheless in the next section we show that in D = 4 the 
simplest 'brute-force' approach works unexpectedly well for the monopole density. 

We also would like to mention that there exists another view on the monopole-artifact 
mixing. Indeed, suppose for a while that the second term in Eq. fll3|) is zero. Then the 
monopole and artifact currents are separately conserved on the hypercubical dual and 
original lattices, respectively, and therefore there is no mixing between them. Thus it is 
the second term in Eq. ([13]) which mixes together the two types of magnetic charges. 

Finally, let us note that there might be another way to fight against the artifacts on 
the lattice. Here we mean in particular a kind of blocking procedure which preserves the 
infrared properties of the theory while filtering out the ultraviolet noise. This might be 
helpful in approaching the limit of infinitesimal lattice spacing where the artifacts for sure 
disappear. 



5. Results of Numerical Simulations. 

In this section we describe the results of our numerical simulations of the four dimensional 
pure SU (2) lattice gauge model. We have noted already that it is in fact crucial to have 



a good approximation to the ground state of the model ([121) which amounts to the ability 
to approximate the global minima of the corresponding action. We have considered three 
different types of minimization procedure, namely, the naive local minimization with and 
without gauge fixing preconditioning and a variant of the simulated annealing algorithm. 

The gauges considered were described above, see Eq. (0), and each gauge was fixed 
by the iterative maximization of the corresponding functional. The gauge was considered 
fixed when all gauge rotation matrices Q x become sufficiently small: 1 — | Tr Q x < 10~ 6 . 
The simulated annealing algorithm which we used is the following. The action (|T3) 
was multiplied by a factor 7 (coupling constant) and then the quantum dynamics of 
variables {z a } was simulated with the standard Metropolis algorithm. From the initial 
value 7 ~ 0.1 the coupling constant was moved towards 7 = 00 with steps A7 = 0.1. 
Fifteen thermalization steps were performed for each 7 and the algorithm was stopped 
when the acceptance rate became smaller than 0.5%. After that on all plaquettes with 
z a = —1 the corresponding i = 0, 3 were inverted n^ a — > —n^ a . 
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To monitor the quality of each algorithm we consider 



n = \d 



(ni, a ■ n jtP ) ) , 



(16) 



where n^ a ■ n Ji( g is the gauge invariant scalar product of two neighboring vectors. The 
dependence r/((3) on 12 4 lattice is presented on Fig. 8, where circles denote the local mini- 
mization without any preconditioning, while the local minimization with minimal Landau 
and Maximal Abelian gauge preconditioning are the squares and diamonds, respectively. 
Triangles are the results for the simulated annealing algorithm. 

It should be noted that although the absolute values of rj varies slightly for various 
algorithms, the corresponding monopole densities differ by an order of magnitude. Note 
also that in the case of Maximal Abelian gauge preconditioning there was in fact no local 
minimization at all (see section 4). Nevertheless the 77 values and monopole densities with 
MAA gauge are at least comparable with ones obtained via the simulated annealing. 

Let us now turn to the physically interesting question about the monopole density. 
Fig. 9 represents the densities of artifacts (circles) and monopoles (squares) versus bare 
coupling constant (3 on a 14 4 lattice. In view of the above discussion of the monopole- 
artifact mixing it should not be a surprise that none of the graphs on Fig. 9 follows the 
renormalization group prediction 
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which is to be valid for physical monopoles. Moreover, Fig. 9 clearly confirms the actionless 
nature of artifacts since their density is almost constant in the whole range of (5 considered. 

As we have already discussed, the monopole-artifact separation is a priori not a well 
defined procedure at a finite lattice spacing. Nevertheless, it is amusing to note that in case 
of D = 4 SU(2) lattice gluodynamics the difference between the monopole and artifact 
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densities follows Eq. flT7|) sharply, see Fig. 10. It seems so that in four dimensions and in 
the weak coupling regime artifacts do indeed form small dipoles with the monopoles living 
on the dual lattice. In this way one can explain the observed scaling behavior of p p h ys — 
Pmonopoies ~ parti facts although there are strictly speaking no rigorous arguments why p phys 
scales. Indeed, the monopole-artifact dipole picture is just an unproved assumption which, 
however, works unexpectedly well in D = 4 SU (2) LGT. The rigorous artifact separation 
is still an unresolved problem which is beyond the scope of the present publication. 

Assuming that the above dipole picture is valid one can estimate the density of the 
physical monopoles in the continuum limit. As follows from Fig. 10, In p p h ys / 4A 3 ~ 11.9 
and using the numerical value of the string tension a^fa = 0.1326 at (3 = 2.6 (see, e.g. 
Ref ||) one gets: 

p phys = (1.72 v^) 3 « (760 MeV) 3 , (18) 
where the conventional value yfa = 440 MeV has been used. 

Conclusions 

We have implemented the gauge invariant monopole charge definition in the SU(2) lattice 
gluodynamics. The lattice implementation is of particular importance since there are no 
analytical approaches to the gauge invariant monopoles in the continuum limit. Indeed, 
the corresponding classical Wu-Yang solution is unstable. Our lattice construction is 
based on the continuum definition of the monopole charge proposed recently and in fact 
follows it as close as possible. In this way we were able to obtain integer valued conserved 
monopole currents. 

As might be expected, the lattice formulation brings in a set of specific problems. 
One of them is the appearance of the effective Ising model. In fact, the gauge invariance 
of monopoles on the lattice crucially depends on our ability to find its ground state, 
which is a non-trivial task because the model is frustrated. But at least in principle this 
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problem might be solved. Moreover, here we found a hint on the distinguished role of the 
Maximal Abelian gauge, which proved to be successful in the context of Abelian monopole 
confinement mechanism. 

Another problem specific for the lattice is the emergence of the lattice artifacts. In 
more detail: 

i) Artifacts are mixed up with physical monopoles. Thus it is necessary to subtract 
their contribution in order to get physical results. For simple quantities, like monopole 
density the artifact separation is indeed possible. Then the density of physical monopoles 
seems to scale correctly towards the continuum limit. However, it is not yet completely 
clear how to get rid of the artifacts in general case. 

ii) Due to artifacts mixing the hypercubical geometry of the original lattice has to be 
changed. The monopole currents are still conserved, but the conservation takes place on 
a complicated lattice, the geometry of which is still not well understood. Probably the 
understanding of the lattice geometry will help to treat artifacts in a more rigorous way. 
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